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Abstract Studying the behaviour of a quantum field in a classical, curved, 
spacetime is an extraordinary task which nobody is able to take on at present 
time. Independently by the fact that such problem is not likely to be solved 
soon, still we possess the instruments to perform exact predictions in special, 
highly symmetric, conditions. Aim of the present contribution is to show how 
it is possible to extract quantitative information about a variety of physical 
phenomena in very general situations by virtue of the so-called Hamilton- 
Jacobi method. In particular, we shall prove the agreement of such semi- 
classical method with exact results of quantum field theoretic calculations. 



1 Introduction 

Suppose we are interested in studying the behaviour of a field ^{x) (scalar, for 
sake of simplicity) in a curved spacetime endowed with a (trapping) horizon 
(e.g. in the vicinity of a black hole). Based on physical intuition, we expect 
that the interaction from the quantum field and the classical background 
gives rise to different phenomena, as: Hawking radiation through the hori- 
zon; decay of unstable particles scattering off the gravitational field; vacuum 
particle creation in regions of strong gravity; radiation from (possibly, naked) 
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singularities, etc. The aformcntioned topics would pertain the investigation 
of a quantum theory of gravity, the lack of which obliges us to work with 
standard techniques. 

The field is governed by the Klein-Gordon equation, 

(□ x - = , (1) 

where the parameter m 2 is interpreted as the field mass and for convenience 
we have inserted ft explicitly. Field quantization can be performed in cither 
canonical or path integral ways: 

<P(x.) — I Dp- Dx • DN x (Gauge-fixing conditions) x exp(z I\p,Z,N]) (2) 

where -A(x) represents appropriate boundary conditions and I[p, x, N] is the 
Hamiltonian formulation of the action. The fact that $(N, x) has to satisfy 
the equation of motion (JTJ imposes constraints on the allowed boundary 
conditions, A. 

Solution to @ is largely unknown due to the difficulty of computing path 
integrals in curved spacetimes. However, some information is accessible in the 
WKB regime of approximation. In this case, one generally finds appropriate 
to look for a solution in the form, 

<P(x) ^ D(x)e" 1 ^ 1 +0(h) (3) 

where the small parameter ft is used to govern the WKB expansion. Inserting 
into ([TJ and equating powers of ft, we obtain to the lowest orders 

ft" 2 : -V a /V a / + m 2 =0 

ft" 1 : 2VD ■ VI + DDI= 

ft : UD = 

.••: ••• (4) 

Exact computation of the pre-factor D(x) is complicated even in very addo- 
mesticated situations and therefore it is beyond our present goal. 
Let us split / into a real and a purely imaginary part: I(x) :— Ir(x) — iS(x); 
then becomes: 

- (V/ K ) 2 + (VS) 2 + m 2 = 0. (5) 

If the imaginary part of I varies with x much more rapidly than the real part, 
that is, if \VS\ ^> |V/r|, it follows from ([S]) that S will be an approximate 
solution to the (Lorentzian) Hamilton- Jacobi equation 



g ab dS a dS b + m 2 = 0. 



(6) 
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Furthmore, the wave function ([3]) will then be predominantly of the form e lS . 
Of course, going from the exact path integral form @ to the approximate 
regime ([3J with / solution to (0}, wc loose specification of the boundary 
conditions A. It will be evident later that this lost is only apparent. 

The basic idea proposed some time ago by Parick & Wilzcek [T] is to in- 
terpret the spacetime horizon - say, for example, of a black hole - as a sort 
of barrier and to study the tunnelling of field quanta through it. Certainly, 
the horizon behavies in quite a different way with respect to usual quantum 
mechanical potential barriers. In ordinary quantum mechanics, the barrier is 
represented by the region between the turning points of the classical trajec- 
tories. Here instead, the horizon is just a point on the classical characteristic 
curves. Evaluating the tunnelling in quantum mechanics means computing 
the ratio between particle wavefunction on the two sides of the barrier. In 
the black hole case, instead, the required ratio is generated by a discontinuity 
in the wavefunction. Moreover, in the familiar context of tunneling through 
a barrier, an imaginary part comes from a negative eigenvalue of the small 
disturbance operator around the classical bounce, while the Euclidean action 
is real. In the black hole case, instead, as we shall see later, it is the action 
itself that is complex. 

Given the whole sort of specifications above, we can conclude the reasoning 
and invoke the well known result according to which, the creation probability 
per unit time of quanta of mass m is given - to leading order in H - by the 
WKB formula 



with S solution to Remarkably, as it has been shown in [3J [3J |U \5\ [5] , 
the procedure outlined so far reproduces the infamous Planckian spectrum 
of Hawking radiation in the case of black hole spacctimes: r oc cxp(— /3ujh), 
with ujh the energy of tunnelling particles through the black hole horizon 
and P interpreted as the inverse temperature of the thermalized field quanta. 
The Hamilton-Jacobi method of tunnelling has therefore proved an elegant 
way to interpret Hawking radiation as a tunnelling process and to derive in 
relatively simple way the associated temperature (T = fi^ 1 ). As we shall 
try to show in the following, the method does not exhaust its power in the 
computation of black hole Hawking temperature, generalizing indeed to a 
wider class of spacetime horizons (e.g. cosmological horizons) and to other 
kinds of semi-classical phenomena. 

We use the conventions according to which the metric signature is (— , +, +, +); 
first latin indices as a, & run over 0, . . . , 3, mid-latin indices as i,j only over 
0, 1. From nown on, we implement natural units, so that c = h = G = fcg = 1- 




(7) 
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2 The Kodama— Hayward formalism for spherically 
symmetric spacetimes 

In the following, we shall limit ourselves to focus only on spherically sym- 
metric spacetimes where no gravitational waves production is involved. The 
line element can be locally written as [7J 

ds 2 = j ij (x)dx l dx j + R 2 (x)df2 2 , (8) 

where the two-dimensional metric jij (x) is referred to as the normal metric, 
{x 1 } are associated coordinates and R(x % ) is the areal radius, considered as 
a scalar field in the two-dimensional normal space. We recall that to have a 
truly dynamical solution, i.e. to avoid Birkhoff 's theorem, the spacetime must 
be filled with matter everywhere. Examples are the Vaidya solution, which 
contains a flux of radiation at infinity, and FRW solutions which contain a 
perfect fluid. 

A dynamical trapping horizon, if it exists, is located at = x( x )\hi with 
x{x) :— 7 4 - 7 (x)diR(x)djR(x), provided that dix(x) = 0. The dynamical sur- 
face gravity associated with the horizon is given by the normal space scalar 
kh = ^\Z\^R(x)\h as proved in [7J. 

In the spherical symmetric dynamical case, it is possible to introduce the 
so-called Kodama vector field K, with {K a G a b)' b = 0, that can be taken 
as its defining property, [8|. It follows that K is a. natural generalization of 
the Killing vector of stationary spacetimes. Given the metric ([8]), the non- 
vanishing Kodama vector components are K l = e l ^djR{x) / (e nl = +1). 
The Kodama vector gives a preferred flow of time and in this sense it gen- 
eralizes the flow of time given by the Killing vector in the static case. As a 
consequence, we may introduce the invariant energy associated with a particle 
of mass m by means of the scalar quantity on the normal space 

lu = -K ■ dS (9) 

where S is the particle action which we assume to satisfy the reduced 
Hamilton-Jacobi equation 

j ij diSdjS + m 2 = . (10) 

Remarkably, the probability rate (Jjj does not depend by the choice of 
coordinates, since the horizon location, the horizon surface gravity, the Ko- 
dama energy are all invariantly defined in the space normal to the spheres of 
symmetry [9]. 

The basic idea can now be roughly described as follows: the reduced 
Hamilton-Jacobi equation (|10[) supplemented by the Kodama energy for- 
mula ([9]) constrains particle's momenta, e.g. d+S — d+S(x :iz , m, u>); thus, the 
mass parameter m gives two complementary energy scales so that, according 
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to the physical phenomenon involved, the two scales exchange the leading 
role in the analysis. 

More in detail, suppose we are interested in the physics of the horizon: tun- 
nelling through the horizon - typically related to Hawking/Unruh effects - 
corresponds to the existence of a simple pole in particle's momenta. In this 
case, it turns out that the mass parameter can be neglected so that, to all 
the extents, particles move along null trajectories. On the other hand, if we 
arc now interested in bulk effects, away from any horizon, then the mass 
parameter plays a crucial role being possibly responsable for a branch point 
singularity in tunnelling particle's momenta. 

Let us make an example in order to make clearer what we mean. As fully de- 
scribed in [TO], the FRW spacetime with spatial curvature k = J| (k = 0, ±1 
and I an opportune length scale) represents a dynamical, spherically symmet- 
ric spacetime exhibiting a cosmological horizon in correspondence of what we 
shall call the Hubble radius, namely Rh (t) :~ (H 2 + fc/a 2 ) -1 / 2 and R{t,r) := 
a(t)r. The Kodama energy is u = yl — kr 2 (—dtJ + rHd r I) = yl — kr 2 uj. 
The Hamilton-Jacobi equation reads —{dtS) 2 + ^ife ^ (d r S) 2 + m 2 = 0, so 
that the radial particle's momentum is 



aHCj(ar) ± a^u 2 - m 2 (l - (ar/flgp) 
1 - {ar/R H y 

Near the horizon, the mass coefficient vanishes so that we can set m = 0. 
Thus, making a null-horizon radial expansion, the action for particles coming 
out of the horizon towards the inner (untrapped) region is S — 2 J drd r I, 
with d r S exhibiting a simple pole at the horizon. To deal with the simple 
pole in the integrand, we implement Feynman's ie-prescription, something 
which resambles the recovering of the boundary conditions encoded in the 
path integral approach mentioned above. Because of ([TJ, r ~ cxp(— ujh/T), 
ujh > for physical particles and T = —kh/2'k {kh < for trapping hori- 
zons of the inner type such as the Hubble radius, cf. [TT] ) the dynamical 
temperature associated to FRW horizon. 

To treat instead the decay of unstable composite particles inside the Hubble 
horizon (i.e., in the untrapped region), we need to identify the energy of the 
particle before the decay as the Kodama energy, w; then we denote by m the 
effective mass parameter of one of the decay products, after the decay. With 
these understandings, we find out that for the unstable particle sitting at rest 
at the origin of the comoving coordinates, there is an imaginary part of the 
action as the decay product tunnels into the region < r < r$ to escape be- 
yond r Q , with r implicitly defined through [a(t)r ] 2 = Rq := (l — R 2 H . 
Assuming a two-particle decay, the rate is 

r = r oe - 27rRH{m - ul} (12) 
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and -To depending on the interaction coupling (e.g. _T ~ A 2 for a Xtfi 3 inter- 
action) . Equation (fT2)l agrees with Volovik result for de Sitter space [H] and 
with asymptotic quantum field theory calculation by Bros et al, |13j . 



3 Vacuum Particle Creation and Emission from Naked 
Singularities 

A perfectly legitimate question we can ask ourselves is weather the method 
is extendable to the case of static black holes as well. With regard to this, 
we consider the exterior region of a spherically symmetric, static, black hole 
spacetime and repeat the same argument. Quite generally, we can write the 
line element as 

rlr 2 

dg 2 = _ e W C{r)dt 2 + ^_ + r 2 d[2 2 (13) 

C(r) 

The analysis of the radial momentum is made easier by setting the Kodama 
energy cj = 0: in the intention, this would correspond to particle creation 
from vacuum, 

r 2 1 

drd r S = m dr (14) 



i-i 



The integration is taken over any interval (ri,r2) where C{r) > 0. Equation 
(fT4|) shows that, under very general conditions, in static black hole space- 
times there could be a decay rate whenever a region where C(r) is positive 
exists. However, it is an easy task to show that the spacelike singularity of the 
Schwarzschild black hole does not emit particles in the semi-classical regime: 
in the interior, the Kodama vector is spacelike, thus no energy can be intro- 
duced there. 

The situation is very different when a naked singularity is present. Consider- 
ing a neutral particle in the Rcissncr-Nordstrom solution with mass M and 
charge Q > (for definitcness) , the line element is 

ds 2 = J r ~ r -^- r +\ t 2 + £ dr ? + r 2 d tf (15) 

r z (r — r_J(r — ?*+) 

with r± functions of (Af, Q) denoting the inner and outer horizons. The 
function C(r) = (r — r_)(r — r + )/r 2 is negative in between the two horizons, 
where the Kodama vector is spacelike, so there the action is real. On the 
other hand, it is positive within the outer communication domain, r > r + , 
but also within the region contained by the inner Cauchy horizon, that is 
< r < r_. Thus, because of (|T4"|) and assuming the particles come created 
in pairs, we obtain that, modulo the pre-factor over which we have nothing to 
say, there is a creation probability per unit time and unit volume (equation 
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P6p not depending upon the creation event) of neutral particles of mass m by 
the strong gravitational field near the Reissner-Nordstrom naked singularity 
(M, Q) which goes as 



At a first look, the process of particle production in the region close to the sin- 
gularity raises the issue of the stability of the solution. However, this does not 
seem to be a problem. In fact, radiation created by the bulk close to the sin- 
gularity comes into the singularity with infinite red-shift and approaches the 
future inner, classically unstable, horizon with infinite blue-shift. Thus, the 
contribution of the radiation coming in the singularity to the back-reaction 
is negligible and the causal structure of the singularity safe; while the blue- 
shifted radiation approaching the future sheet of the inner horizon will con- 
tribute to its quantum instability (Cf. [T3] for further investigation). 

A complementary and potentially interesting effect is the emission from 
the naked singularity itself. We investigate this problem for the case of two- 
dimensional dilaton gravity, and will come back to Reissner-Nordstrom solu- 
tion afterward. 

Consider the two-dimensional metric 

ds 2 = o~ x dx + dx~ , a := \ 2 x + x~ — a(x + — Xq)9(x + — Xq) (17) 

where A is related to the cosmological constant by A = — 4A 2 and a represents 
the wave amplitude. This metric arises as a solution of 2D dilaton gravity 
coupled to a bosonic field with stress tensor T ++ = 2a S(x + — Xq ), describing 
a shock wave, a = is a naked singularity partly to the future of a flat space 
region (linear dilaton vacuum). The heavy arrow in the figure represents 
the history of the shock wave responsible for the existence of the timelikc 
singularity. The goal is to compute the outgoing flux which is given by: 2T ++ — 




(16) 



wave. 



Fig. 1 The naked singu- 
larity formed by the shock 
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2T In order to do this, we notice that the Hamilton- Jacobi equation 

implies either d+S = or d^S = 0, S being the action. To find the ingoing 
flux we integrate along x + till we encounter the naked singularity, using 
d-S = 0, so that S = J dx + d+I = J dx + — . The absorption probability 
is r(u>) = i~b cxp[— 2itlu/(\' 2 x~ — a)}. The flux is computed by integrating 
the probability over the coordinate frequency u> = ^ (that is, the variable 
conjugate to the time coordinate), with the density of states measure 

rp (\ 2 x~— a) 2 

1 ++ - 1 2 ^a'\ ■ 

To find T , we integrate now along x starting from the naked singularity, 

this time using d+S = 0. A similar calculation gives T = Jq \l*i a l ■ 

is given by the conformal anomaly, T = 4aT^ = R/24ir (for one bosonic 

degree of freedom). Matching to the anomaly gives To = 7r 2 /24 ~ 0(1). 
Note that the stress tensor diverges approaching the singularity, indicating 
that its resolution will not be possible within classical gravity but requires 
instead quantum gravity [15]. 

Indeed, all these results agree with the one-loop calculation to be found in 

Returning now to the Reissncr-Nordstrom solution, could it be that the 
naked singularity emits particles? In the four-dimensional case one easily sees 
that the action has no imaginary part along null trajectories either ending 
or beginning at the singularity. Formally this is because the Kodama energy 
coincides with the Killing energy in such a static manifold and there is no 
infinite red-shift from the singularity to infinity. Even considering the metric 
as a genuinely two-dimensional, however, it is possible to show that the action 
does not exhibit any imaginary part |10j . 

It is fair to say that the Reissner-Nordstrom naked singularity will not emit 
particles in this approximation something which seems to be coherent with 
quantum field theory results, [TU1 IT7]. 



4 Conclusions 

We have shown that semi-classical tunnelling method can handle several 
quantum effects: radiation from dynamical horizons (both cosmological and 
collapsing); gravitational enhancement of particle decay otherwise forbidden 
by conservation laws; radiation from two-dimensional naked singularities. 
Normally, great efforts are needed to analyze quantum effects in gravity, while 
instead the tunneling picture promptly gives strong indications of what's go- 
ing on. The obtained agreement between both the particle decay rates and the 
radiation from naked singularities in the tunnelling picture and the (asymp- 
totic of the) exact results - when they exist in particular conditions - gives 
confidence, in our opinion, of the validity of the method even in more general 
situations. 
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